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1 Introduction 



In this paper we consider pseudo-Euclidean Toda-like system described by the following 
Lagrangian 

L = L{z\z\M) = ^tf- 1 TjabZ a Z b -jW(z), (1.1) 

where M > is the Lagrange multiplier (modified lapse function), 

{Vab) = diag(—l, +1, . . . , +1) is matrix of minisuperspace metric, a, b — 0, . . . , n — 1, and 

rn 

V» = £A Q exp«z a ) (1.2) 

a=l 

is the potential. Here A a ^ 0. 

Models of such type occur in multidimensional cosmology (see, for example, PQ-[H] and 
references therein) as well as in multicomponent 4-dimensional cosmology (Tj-pJJ]. The 



Lagrange systems (1.1), (1.2) are not well studied yet. We note, that the Euclidean Toda- 
like systems are more or less well studied [12]- [18] (at least for certain sets of parameters, 
associated with finite-dimensional Lie algebras or affine Lie algebras). There is also a criterion 
of integrability by quadrature (algebraic integrability) for the Euclidean Toda-like systems 
established by Adler and van Moerbeke [16]. We note that in gravitational context the 



Euclidean Toda-like systems were first considered in [T9]-|21j. 

In this paper we consider the behavior of the dynamical system (1.1) for n > 3 in the 
limit 

z 2 = -(z ) 2 + (z) 2 — > — oo, z = (z°, z) G V_, (1.3) 
where V_ = {(z°, z) G R n \z° < —\z\} is the lower light cone. The limit (1.3) implies 

z° -> -oo (1.4) 

and under certain additional assumptions describes the approaching to the singularity 
in corresponding cosmological models. We impose the following restrictions on the vectors 



in the potential (1.2) 



1) A a > 0, if {u a f = -K) 2 + {u a f > 0; (1.5) 

2) u% > for all a = l,...,m. (1.6) 

We note that in multidimensional cosmology a special interest is connected with the 
investigations of oscillatory behavior of scale factors near the singularity [22]-[26]. This 
direction in higher-dimensional gravity was stimulated by well-known results for "mixmaster" 
model [7]- [10]. We note, that there is also an elegant explanation for oscillatory behavior of 
scale factors of Bianchi-IX model suggested by Chitre [SI 110] and recently considered in [2"7] - 
|29j . In the Chitre's approach the Bianchi-IX cosmology near the singularity is reduced to a 
billiard on the Lobachevsky space H 2 (see Fig. 4 below). The volume of this billiard is finite. 
This fact together with the well-known behavior (exponential divergences) of geodesies on 
the spaces of negative curvature may lead to a stochastic behavior of the dynamical system 
in the considered regime [301 EI] ■ 
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Here we consider the generalization of Chitre's approach [S] to the multidimensional case 
|26j . In the limit (1.3) the dynamics of the model (1.1 ), ( 1.2[ ) (with the restrictions (1.5), (1.6) 
imposed) is reduced to a billiard on the (n — l)-dimensional Lobachevsky space H n ~ l (Sec. 
2). The geometrical criterion for the fmiteness of the billiard volume and its compactness is 
suggested. This criterion reduces the considered problem to the geometrical (or topological) 
problem of illumination of (n — 2)-dimensional unit sphere S n ~ 2 by m + < m point-like 
sources located outside the sphere [32} [33]. These sources correspond to the components 
with (u a ) 2 > (Sec. 2). When these sources illuminate the sphere then, and only then, 
the billiard has a finite volume and the corresponding cosmological model possesses an 
oscillatory behavior near the singularity. In this case from topological requirements we 
obtain the restriction on the number of components with (u^) 2 > 0: m + > n, i.e. this 
number should be no less than the minisuperspace dimension. In Sec. 3 we illustrate the 
suggested approach using the examples of Bianchi-IX and prototype cosmological models. 



It should be noted that fixing the gauge in (1.1 ) Af = 1 we get Lagrange system 

1 



Li 



Vab z a z b -V(z), 



with the zero-energy constraint 



Ei 



~-r] ab z a z b + V{z) = 0. 



For non-zero energies E\ we are lead to Lagrange systems with the Lagrangians 



1 

2 L 



L x + ~-e{z n ^ 



;i.t) 



x\ 



;i.9) 



—signal, and the zero-energy constraint E 2 = 0. Thus, for E\ < we get the 



pseudo-Euclidean Toda-like system (1.1) with n replaced by n + 1 and the potential (1.2). In 



cosmology this corresponds to addition of scalar field. For E\ > we obtain the Lagrangian 
system with the kinetic term of signature (—,+,...,+,—). In cosmology this corresponds 
to the incorporation of a phantom field. We also note that the Euclidean Toda-like systems 
may be embedded into (1.1) by considering the potentials (1.2) with Uq = 0. 



2 Billiard representation 



Here we consider the behavior of the dynamical system, described by the Lagrangian (1.1) 



for n > 3 in the limit (1.3). We restrict the Lagrange system (1.1) on V_, i.e. we consider 
the Lagrangian 

L_=L| rM _, M_ = V_xi? + , (2.1) 

where TM_ is tangent vector bundle over M_ and R + = {Af > 0}. (Here F\a means the 
restriction of function F on A.) Introducing an analogue of the Misner-Chitre coordinates 

in v_ p rug 

i + y 2 
, y 



-exp{-y 



-2 exp(— y 



i-y 



TV 



(2.2) 
(2.3) 
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\y\ < 1, we get for the Lagrangian (1.1] 



IaTV^H?/ ) 2 + h l3 {y)ijY]-UV. 



Here 

i,j — 1, . . . , n — 1, and 
where 



2\-2 



a=l 



<%, u) = -e-y (1 - fy'luoil + + 2uy\, 
We note that the (n — l)-dimensional open unit disk (ball) 

D n - 1 = {y = (y\...,y n - 1 )\\y\<l}cR n - 1 



(2.4) 
(2.5) 

(2.6) 
(2.7) 

(2.8) 



with the metric h = hij(y)dy l Cg> dyi is one of the realization of the (n — 1) -dimensional 
Lobachevsky space H n ~ x . 
We fix the gauge 

TV = exp(-2y°) = -z 2 . (2.9) 

Then, it is not difficult to verify that the Lagrange equations for the Lagrangian ( 1.1[ ) with 
the gauge fixing (2.9) are equivalent to the Lagrange equations for the Lagrangian 



with the energy constraint imposed 



(2.10) 

(2.11) 

(2.12) 
(2.13) 

Now we are interested in the behavior of the dynamical system in the limit y° — > — oo 



Here 



where 



E* = - l -(y ) 2 + \h l3 {y)y l V 3 +V m = 0. 



in 

V* = e~ 2y °V = ^«exp($(?/,M Q )), 

a=l 



${y,u) = -2y° + <$>{y,u). 



(or, equivalently, in the limit (1.3) implying (1.4). Using the relations (u ^ 

•fc.)=-«.« f <-V)^P-V, 



we get 



A(y,v) = (y-v) 2 -v 2 + l 



lim exp $(?/, u) = 

y —*— oo 



(2.14) 
(2.15) 

(2.16) 



for u 2 = — ul + (u) 2 < 0, u > and 



for u 2 > 0, Mo > 0. In (2.17) we denote 



lim exp $(y, tt) = ^(-^(f, -u/u )) 



(2.17) 



^oo(x) = + OO, X > 0, 

0, x < 0. 



Using restrictions (1.5), (1.6) and relations (2.12), (2.16), (2.17) we obtain 



lim V,(y°,y) = E ^oo(-A(y,- M «K))- 



y — * — oo 



oeA 4 



Here we denote 



A + = {a|K) 2 >0}. 
The potential Voo may be also written as following 

V oo (0) = V(y,B)= 0, yEB, 

+oo, yeD^XB, 

where 

B = f| 5(m q ) C D"- 1 , 

aeA + 



B{u a ) = I y E D 



n-l 



u 



> 



\ 



It 



(2.18) 

(2.19) 
(2.20) 



(2.21) 



(2.22) 



(2.23) 



a E A + . B is an open domain. Its boundary OB = B \ B is formed by certain parts of 
m + = |A+| (m + is the number of elements in A + ) of {n — 2)-dimensional spheres with the 
centers in the points 



v° = -tT/uJ, a E A + , 



(|w a | > 1) and radii 



respectively (for n — 3, m + = 1, see Fig. 1). 

So, in the limit y° —>■ —oo we are led to the dynamical system 



L oo = -\{y ) 2 + \h ] {y)y^ -y K 



V) 2 + fM£h/V + Kx 



(2.24) 



(2.25) 



(2.26) 
(2.27) 



which after the separating of y° variable 
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y° = u(t-to), (2.28) 
(a? 7^ , to are constants) is reduced to the Lagrange system with the Lagrangian 

1 



JB = My)y l y 3 -v^B). 



Due to (2.28) 



E f 



-h i3 {y)ijY +V(y,B) 



LO 

T 



(2.29) 



(2.30) 



We put L) > 0, then the limit t — 
empty (A+ = 0), we have B = D 71 ^ 1 and the Lagrangian (2.29) describes the geodesic flow 
on the Lobachevsky space if n_1 = (D™ _1 , lively 1 ® 
of non-trivial geodesic solutions (i.e. y(t) ^ const): 



oo corresponds to (1.3). When the set (2.20) is 
mgian (2.29) describes the geodesic flow 
dyi). In this case there are two families 



1. y(t) = n\[y/v 2 — 1 cos ip(t) — v] + n 2 \/v 2 — 1 sin (p(t), 

(p(t) = 2 arctan[(t> — \Jv 2 — 1) tanh(|o;(t — ti)], 

2. y(t) = n 2 tanh(|w(t-ti)). 



(2.31) 
(2.32) 
(2.33) 



Here n\ = n\ = 1, n\n-i — 0, v > 1, a> > 0, ti = const. 

Graphically the first solution corresponds to the arc of the circle with the center at point 
(— vfii) and the radius \/v 2 — 1. This circle belongs to the plane spanned by vectors n\ and 
f?2 (the centers of the circle and the ball D n ~ l also belong to this plane). We note, that the 



solution (2.31 )-(2.32) in the limit v — * oo coincides with the solution (2.33). 

When A + ^ the Lagrangian (2.33) describes the motion of the particle of unit mass, 



moving in the (n — l)-dimensional generalized billiard B C D n (see (2.22). The geodesic 



motion in B (2.31 )-(2.33) corresponds to a "Kasner epoch" and the reflection from the 
boundary corresponds to the change of Kasner epochs. For n = 3 some examples of (2- 
dimensional) billiards are depicted in Figs. 2-4. 
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Figure 2: Billiard with infinite volume for n — 3, m + = 3. 




Figure 4: Billiard corresponding to Bianchi-IX model (non-compact with finite volume). 



We note, that the boundary of the billiard dB is formed by geodesies. For some billiards 
this fact may be used for "gluing" certain parts of boundaries. 

The billiard B in Fig. 2. has an infinite volume: volB = +00. In this case there are 
three open zones at the infinite circle \y\ = 1. After a finite number of reflections from 
the boundary the particle (in general position) moves toward one of these open zones. For 
corresponding cosmological model we get the " Kasner-like" behavior in the limit t — > —00 

For billiards depicted in Figs. 3 and 4 we have volB < +00. In the first case (Fig. 
3) the closure of the billiard B (in the topology of D n_1 ) is compact (in the topology of 
D"^ 1 ) and in the second case (Fig. 4) B is non-compact. In these two cases we have an 
"oscillatory-like" motion of the particle. 

Analogous arguments may be applied to the case n > 3. So, we are interested in the 
configurations with finite volume of B. 

We propose a simple geometric criterion for the finiteness of the volume of B and compact- 
ness of B in terms of the positions of the points (2.24) with respect to the (n— 2)-dimensional 
unit sphere S n ~ 2 (n > 3). 

Proposition 1. The billiard B (2.22) has a finite volume if and only if the point-like 
sources of light located at the points v° (2.24) illuminate the unit sphere S n ~ 2 . The closure 
of the billiard B is compact (in the topology of D"^ 1 ~ H n ~ r ) if and only if the sources at 
points (2.24) strongly illuminate S n ~ 2 . 

Definition. Here the point is called strongly illuminated if it belongs to the interior 
(open) part of the illuminated region. At Fig. 1 the source P strongly illuminates the open 
arc (Pi,P 2 ). 

Proof. We consider the set d c B = B c \ B, where B c is the completion of B (or, 
equivalently, the closure of B in the topology of R n ~ l ). We recall that B is the closure of B 
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in the topology of D n 1 . Clearly, that d c B is a closed subset of S n 2 , consisting of all those 



points that are not strongly illuminated by sources (2.24). There are three possibilities: i 



d c B is empty; ii) d c B contains some interior point (i.e. the point belonging to d c B with some 
open neighborhood); iii) d c B is non-empty finite set, i.e. d c B = . . . y{\. The first case i) 
takes place if and only if B is compact in the topology of D 11 ^ 1 . Only in this case the sphere 
S n ~ 2 is strongly illuminated by the sources (2.24). Thus the second part of proposition is 



proved. In the case i) volB is finite. For the volume we have 



volB= [ d n ~ 1 yVh= [ dr(l -r 2 f~ n S T . 

JB JO 



(2.34) 



The "area" S r — > C > as r — > 1 in the case ii) and, hence, the integral (2.34) is 



divergent. In the case iii) 



S r ~ Ci(l - r 



2(n-2) 



as r 



1 



(2.35) 

(Ci > 0) and, so, the integral (2.34) is convergent. Indeed, in the case iii), when r — > 1, 
the "area" S r is the sum of / terms. Each of these terms is the (n — 2)-dimensional "area" 
of a transverse side of a deformed pyramid with a top at some point y^, k = 1, . . . , / . This 
multidimensional pyramid is formed by certain parts of spheres orthogonal to S n ~ 2 in the 
point of their intersection y^. Hence, all lengths of the transverse section r = const of the 
"pyramid" behaves like (1 — r) 2 , when r — ► 1, that justifies (2.35). But the unit sphere S n ~ 2 
is illuminated by the sources (2.24) only in the cases i) and iii). This completes the proof. 



To our knowledge, the problem of illumination of convex body in a vector space by point- 
like sources for the first time was considered in [321 E3]- For the case of S n ~ 2 this problem 
is equivalent to the problem of covering the spheres with spheres (3U [35]. There exists a 
topological bound on the number of point-like sources m + illuminating the sphere S n ~ 2 [33J: 



m + > n. 

Thus, we are led to the following proposition. 

Proposition 2. When m + < n, i.e. the number of components with 



(2.36) 



V 1 



) 2 > is less 

than the minisuperspace dimension, the billiard B (12.221) has infinite volume: volB = +oo. 



In this case there exist an open zones on the sphere S n 2 and the oscillatory behaviour 
near the singularity for cosmological models is absent (we get a Kasner-like behaviour for 
t -> -oo 01 ESI). 

Remark 1. L et the points (2.24) form an open convex polyhedron P C R n ~ 1 . Then the 
sources at (12.241) illuminate S n ~ 2 , if D n ~ x C P, and strongly illuminate S n ~ 2 , if D n ~ l C P. 



3 Some examples 

3.1 Bianchi-IX cosmology 

Here we consider the well-known "mixmaster" model [7J |S] with the metric 



exp[27(t)]dt dt + J2exp[2x i (t)]e i ® e { 



(3.1) 
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where 1-forms e % = e l u (()d( u satisfy the relations 

de l = -e ijk e j A e k , (3.2) 

i,j, k = 1, 2, 3. The forms e % are defined on S 3 ~ SU(2) and are the components of the 
Morera-Cartan form (see, for example [36J) on SU(2). The Einstein equations for the metric 



(3.1) are equivalent to the Lagrange equations for the Lagrangian 



L = L( 7 , x, x) = ) 2 U- l G ij x i x j - MV{x). (3.3) 



Here 



Gij = Sij - 1 (3.4) 
are the components of the minisuperspace metric, 

V = V mix = ^(e 4 * 1 + e 4x2 + e 4x3 - 2e 2xl+2x " - 2e 2x2+2x3 - 2e 2xl+2x ") (3.5) 
is the potential and the Lagrange multiplier is 

M = exp( 7 - 7o (x)) > 0, (3.6) 

where 

n 



to =xy ( 3 - 7 ) 



1=1 



and n — 3. 

In z- variables [1,2] 



z o = q -iJ2x\ q=[n/(n-l)} 1 / 2 , (3.8) 

i=l 

z a = [i/(„ - a + l)(n - a)} 1/2 J2 ~ 0> ( 3 - 9 ) 

j=a+l 



a = I, ... ,n — 1, where here n = 3, we get the Lagrangian (1.1 ) with 3- vectors 



u 



1 . . ^ o 1 . . ,.. , 4 



V l>l,-V3), « 2 = -^(l,l,+V3), m 3 = —(1,-2,0), (3.10) 

u 4 = hu 1 + it 2 ) , = -fa 1 + m 3 ), = ^ + A (3- 11 ) 
and, consequently, 

K) 2 = 8, ( M 3+a ) 2 = 0, (3.12) 
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a 



1,2,3. Thus, the conditions (1.5), (1.6) are satisfied. The components with a 



4,5,6 do not "survive" in the approaching to the singularity (see (2.16)). For the vectors 



(2.24) we have 



v 1 



;i,-V3), 



;i,W3), 



-2,0), 



(3.13) 



(see also [27] )■ In this case the circle S 1 is illuminated by 
= 1,2,3, but is not strongly illuminated. In agreement with 



i.e. a triangle from Fig. 4 
sources located at points v 1 , 1 
Proposition 1 the billiard B has a finite volume, but B is not compact. 



3.2 Prototype model 



Now, we consider the model, described by Lagrangian (3.3) with minisupermetric (3.4), 
L, . . . n, and the potential 



1 n 
V = - ]T exp[2 ]T x l + 2(x i - x j - x k )} 

(i,j,k)eS i=l 



(3.14) 



where S = {(i,j,k)\i,j,k = l,...,n;j ^ i ^ k, j < k}, n > 2. This model may be 
considered as a prototype cosmological model describing the behaviour of the solutions to 
the Einstein equations in the dimension D = 1 + n near the singularity [221 123] . The sum in 
contains n(n — l)(n — 2)/2 terms. In ^-coordinates (3.8), (3.9) we get the Lagrangian 



with 



u 



a\2 



8 > and Mq = 2/q > for any component a G S. Thus, the restrictions 



and (1.6) are satisfied. 



The corresponding billiard is not compact for all dimensions, and has a finite volume 
for n < 10 and infinite volume for n > 10 [38J. This proposition follows from Proposition 1 
formulated in terms of inequalities on Kasner parameters [38] and analysis of this inequalities 
performed in Ref. [23]. For n = 4 the billiard was studied recently in [37] . 



3.3 Scalar field generalization 



Let us consider the Lagrangian (J 1 . 1 [) with a, b 
not depend upon the 



0, . . . , n and u% = 0, i.e. the potential does 



"scalar field" ip = z n . In this case the correponding billiard has an 
infinite volume, since at least two points (0, . . . , 0, ±1) ("North and South Poles") on the 
"Kasner sphere" S n ~ l are not illuminated by the sources at the points v a with = 0. Thus 
the oscillatory behaviour in this case is absent. We may expect that the points belonging to 
the domain of shadow on the sphere S 1 ™ -1 are integrals of motion in this case. Due to (1.9) 



this analysis may be also applied to pseudo-Euclidean Lagrange systems (1.7) in the region 
of negative energy Ex < 0. 



4 Discussions 



Thus, we obtained the "billiard representation" for pseudo- Euclidean Toda-like system (1.1 ), 



(1.2) and proved the geometrical criterion for the finiteness of the billiard volume and the 
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compactness of the billiard (Proposition 1, Sec. 2). This criterion may be used as a rather 
effective (and universal) tool for the selection of the cosmological models with an oscillatory 
behavior near the singularity. 

Remark 2. It may be shown that the condition (1.5) may be weakened by the following 

one 

< > 0, if {u a f < 0. (4.1) 



In this case there exists a certain generalization of the set B(u a ) from (2.23) for arbitrary 
u o (X^") 2 > 0)- The Proposition 1 (Sec. 2) should be modified by including into consideration 
the sources at infinity (for Uq = 0) and "anti-sources" (for Uq < 0). For "anti-source" the 
shadowed domain coincides with the illuminated domain for the usual source (with Uq > 0). 
In this case we deal with the kinematics of tachyons (v 01 are the velocities of tachions). 
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